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^ Abstract 
oo 

■ The effect of the inevitable couphng to external degrees of freedom of a quan- 

■ turn computer are examined. It is found that for quantum calculations (in 
which the maintenance of coherence over a large number of states is impor- 

• tant), not only must the coupling be small but the time taken in the quantum 

,1^ • calculation must be less than the thermal time scale, h/ksT. For longer times 



X 



the condition on the strength of the coupling to the external world becomes 



. much more stringent. 



I. INTRODUCTION 

Quantum computers have recently raised a lot of interest. A number of papers have 
argued that quantum computers can solve certain problems much more efficiently than can 
classical computers. Shor ||^ has shown that a quantum computer could solve the problem 
of finding discrete logs (mod N) and of finding the factors of a large number N in a time 
which is a polynomial function of the length L of the number. For factoring the best 
known algorithm, the Number Field Sieve takes a time of order exp(c(L)^/'^(ln(L))^/^, 
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where c{L) is roughly constant and approximately equal to 2 for large L. Although this is 
subexponential, it is worse than any polynomial for large N. A crucial feature of the ability of 
quantum computers to be more efficient in certain problems involves having the computer be 
placed in the coherent superposition of a very large number (exponential in L) of "classical 
states" , and having the outputs interfere in such a way that there is a very high probability 
that on the appropriate reading of the output, one would obtain the required answer. One is 
replacing exponentiallity in time with exponentiallity in quantum coherence. This requires 
that the computer be able to maintain the coherence during the course of the calculation. 
This paper examines this requirement, and examines the constraints placed on the ability to 
maintain this coherence in the face of coupling to external heat baths. Landauer |^ has long 
emphasized the necessity of examining the effect of both imperfections and of the coupling 
to the external world of any realistic device on the ability of quantum computers to realize 
their promise. This paper is thus a first step in that direction. 



II. DECOHERING NOISE 

I will look at only the simplest model, in which I ask about the maintenance of coherence 
in a memory of length L. This does not take into account the effect that the course of the 
computation itself would have on the rate of loss of coherence, but I would expect that only 
to increase the problem. Thus let us assume that that the number is represented in the 
computer as a string of binary digits of length of the order of L = ln{N). The memory 
cells will each be taken to be two level systems, with each of the two levels having the same 
energy. The two states will be take to be the eigenstates of a "spin" operator a^. 

In a conventional computer, the way in which the calculation is "kept on track" is by 
including dissipation in order to damp out any attempt by the system to make a transition 
( except of course those driven by the computation) [Q. I will therefore assume that the 
interaction with the environment has the two desired eigenstates of the memory as eigen- 
states of the interaction. The environment will be modeled by a massless scalar field |^ 
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derivatively coupled to the memory cell, so that the the full Hamiltonian is 

H^l ((7r(x) + eh{x)a,f + {d,cf>{x)f) dx (1) 

(The associated lagrangian has the simple derivative coupling form 

L^ll {{dtcl>f - {d^f - 2eh{x)cj>{x)a,) dx.) (2) 

Here h{x) is some interaction range function, and tt is the momentum conjugate to 
The Heisenberg equations of motion are 

TT = (3) 

<j) = n + eh{x)a^ (4) 

The exact solutions for the Heisenberg equations of motion for are 

<P{t, x) = 1 (0(0, x-t) + 0(0, x + t)+ 7r(0, y)dy) (5) 
- i / Wzit y\) - cr{-\x -y- t\)e{x -y-t) 

-(T,{-\x -y + t\)Q{-{x -y + t))] h{y)dy 

where Q{x) = {0 i/ x < 0; 1 if x > 0}. 

Since in the model, cx^ is a constant of the motion (recall that 1 am not taking into 
account the effects of the operation of the computer) the solution for is thus 

0(i, x) = i (0(0, x-t) + 0(0, x + t) + 7r(0, y)dy) (6) 

- la, J[l-Q{x-y-t)- Q{-{x -y + t))] h{y)dyh{y)dy (7) 

I will however be working in the Schroedinger representation in the following. 

L assume that the initial state of the environment is a thermal density matrix Rt with 
temperature T, and the initial state of the spin is a density matrix p(0). The total state 
is assumed to be a product state of these two initial states. The reduced state of the spin 
system at any time (t) after tracing out over the state of the environment is a density matrix 
given by 
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p(t) = l(l + p(t)-a) (8) 

where p{t) is a time dependent vector of length less that or equal to unity. p{t) is given by 

p{t) = Tr (ae'^^p{0)RTe-'^*) (9) 

where the trace is over all of the degrees of freedom of spin system and bath. 
We can write H as 

^ _ g« J eh{x)(j>{0,x)dx(72 jj^^-i J eh{x)(l){0,x)dxaz {W) 

since J h{^)<l>{o,x)dx jg j^^^ translation operator taking 7r(0, a;)to 7:{0,x) + 

e J h{x)(l){0,x)dxaz, and since az commutes with Hq. 
Thus 

p{t) = Tr (^ffe'l't'^^^'^'^'"'^e-'f^^^'^''"^'"''p{0) ■ a (11) 

Xe^^°^ Rxe~^^°^e,^ /0(t)e/i(ia;a-2g-i J 4>{0)ehdxaz^ 

where = e'''^°*0(O, x)e~'^^°^ is the time development of the free field with the same initial 
conditions 0(0) and 7r(0), i.e., 

^(^t,x) ^ -{(f){0,x-t) + (l){0,x + t)+ 7r{0,x')dx') (12) 
2 Jx-t 

Using (T^ — 1 and the fact that Rt is diagonal in the energy representation, we can write 

p{t) as 

p{t) = Tr (^^e^/('^(°)-'^(*»^'^'^^"^ p(0) ■ <ji?ye*/W*)-'^(o))^'^<^^^-) (13) 

(Note that the extra terms from the Cambell-Baker-Hausdorf formula cancel out.) This 
can furthermore be written as 

p{t) = Tr (aim • (^ - (1 - cos{J{4>{t) - m)^hdx))a,e, (14) 
+sin{J{4>{t) - (l){0))ehdx)% x aRr) 



where is the unit vector in the z direction. Because Rt is symmetric in and tt, the sin 
term is zero, and 

J{t) = Tr{RTCOs{ J {4>{t) - (l){0))ehdx))) = e''^ J Tr{Rr{m-m))ehdx (^5) 

We are thus left with 

p,{t) = pM (16) 
Px{t) = J{t)pM (17) 

Py{t) = J{t)Py{0) (18) 

For later use, let us examine J{t) in various regimes. Let us take h{x) such that h{k), the 
Fourier transform of h{x) is of the form e"^^'^. F is a cutoff parameter typical of interactions 
with the environment. I will assume that F >> 1/T. We then get 



ln{J{t)) 



4 / " cos{kt))e-^^^ dk (19) 



We can approximate coth{x) ^ 1 + e ^(^ + ...). This gives us 

Hm) « (\ln (£1±^) - hn (1 + (2Ttr) - mm (i^)) (20) 

There are essentially three regimes for the time dependence of J{t) given by the conditions 
t < r, r < t < 1/T and t > 1/T. In the first regime, t < F, we have approximately 

^^(■^W ^ ^ (21) 
For the intermediate regime, F < t < 1/T ,the quantum regime, we have 

HAt)) ^ -^In (^) (22) 

Finally, for the long time regime t » 1/T, the thermal regime, we have 

ln(J(t)) ^ -e^Tt (23) 

The important feature of these asymptotic formula is that for the intermediate regime, 
which I call the quantum regime since the behaviour is dominated by the vacuum state of 



the environment, In (J) increases only logarithmically with t. In contrast, the third regime, 
the thermal regime, it increases linearly with t. This will be important in determining the 
ultimate size of a number which can be say factored with a quantum computer, because of 
the dependence of the computing time on the length of the number being factored. 

This was for the most familiar case of an "ohmic" coupling to the heat bath. In the case 
of superohmic {h{k{(jj)) = cu^e^ — Tu)) for s > 0), the function ln{J{t)) is essentially constant 
for times less than 1/T and grows as t^~^ in the thermal regime for s < 1. For s > 1, J is 
constant in both regimes, although it is smaller in the thermal regime than in the quantum 
regime. ( and is essentially constant even for such times if s > 1) In the subohmic case, 
— 1 < s < 0, on the other hand, lnJ{t) grows roughly as in the quantum regime and 
as t^~^ in the thermal regime. Again, in the thermal regime the growth in decoherence is a 
factor of t larger than in the quantum regime. 

The above analysis was carried out for a single bit in the memory of the quantum 
computer. Let us examine the situation in which our memory has some large number L of 
bits. Each bit is assumed to couple to its own heat bath of exactly the above type. The 
question now is "What is the rate of of loss of coherence of a coherent sum of numbers stored 
in the memory". le, define the state \n >— > |ni_2 > ■■■\no >, where rii is the ith bit 

of n. Consider a coherent state 

>^^an\n> (24) 

n 

The probability that after time t the memory remains in the the state ip is given by 

Prob^ = < ^\Trenvironment (e'S^^^*|V' > |0 >< 0| < jjj\e'^i'''') (25) 

= ^ a*^an'ama*^,Y[Trenvironi (< ni|e*^'Vi >< ^'ile'^^'^H >) 

nn'mm' i 
nn' i 

where (n, (8) n'^) is the XOR of the ith bits of n and n'. 

This expression tells us how the coherent sum over the various states of the memory 
representing various numbers decoheres as a function of time. As an example, let us chose the 
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completely coherent state in which each of the numbers of length L has an equal probability. 
This state is typical of the state required in performing quantum calculations of the sort in 
which a quantum computer is much faster than a classical computer. le, I choose IttnP = 2~^. 
Furthermore let me assume that each bit is coupled to the environment in exactly the same 
way so that Ji{t) — J{t). Then we have the probability that the coherence will be maintained 
over time t as 

Proh = 2-2^ n Jitf"'®"''^- (26) 

nn' i 

To evaluate this first fix the number n. The number of numbers n' which differ from n in 1 
bit is L. The number which differ in 2 bits is L[L — l)/2 and the number which differ in r 
bits is -uT^,. Thus the above becomes 

Prc^ = E E = 2- Ed + Jf ^ (^)' (27) 

If we assume that 1 — J is very small (which is the only case in which the system has 
any hope at all of acting like a quantum computer) , this is well approximated by 

Proh fti e-^^^^--^^ (28) 

as long as L(l - J) < 1/(1 - J). 

The strength of the quantum computer is that the time required to perform the calcu- 
lation is a polynomial in the length L of the number. This time I will designate by t{L). 
Since the quantum calculation is polynomial in L we can write t(L) a; for a > 1. We 
thus have that the probability of maintaining coherence over the time of the calculation is 
of the order of 

hi{Proh) ^ -0(l)Le^ ln(r(iV)) ^ -0{l)e^L ln(L) (29) 
in the quantum regime while it is of order 

\n{Prob) ^ -0(l)L"+^e2 (30) 
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with a smooth transition between the two regimes. In order to have a reasonable probabihty 
of obtaining the correct answer, one needs the probabihty of obtaining the quantum coherent 
answer to be of order 1. This imphes that one must have a sufficiently small e^, the coupling 
parameter between the heat bath and the system. As long as one is in the quantum regime, 
the relation between the coupling and the maximum length of the number one can handle 
is essentially inverse linear, no matter what the polynomial dependence of the calculation. 
However, once one has entered the thermal regime, a decrease in the coupling buys one only 
a small increase in the length of the number L that one can use. I.e., in the presence of a 
coupling to the heat bath, the thermal time scale ^ = plays a crucial role. As long as 
the calculation can be completed in a time less than this, one can imagine decreasing the 
coupling to the heat bath for the memory cells so as to achieve the maximum L. If however 
the time for the calculation is longer than the thermal time scale, it becomes very difficult 
to decrease the coupling to the bath sufficiently to achieve the necessary coherence. 

Is it possible to use the computer even if the quantum state looses coherence? I cannot 
answer this in general, but can show that one strategy does not work. One could imagine 
trying to make up for the loss of coherence by increasing the number of times the program 
is run. (This is in fact a crucial factor in the Shor algorithm for factoring, not because 
of decoherence, but because the calculation itself has a finite probability of not giving the 
correct outcome.) After a sufficient number of attempts, one should by chance have a system 
which has maintained coherence. In the factoring problem, one can test ones answer ( does 
it give the factors of the number), and simply keep repeating the experiment until one gets 
the right answer. However, in M trials, the probability of never finding a coherent outcome 
is (1 — Proh)^ ^ The number of trials required to make this small (i.e., so that 

one has a high probability of having had a coherent run) is thus , the required number of 
attempts is M fa 1/Prob ~ e!^i'^)^^'^i^) in the quantum regime, which is exponential in the 
length. In the thermal regime, this time scale is even worse. One will thus have lost all 
advantages of the quantum nature of the computer. We see that one must make sure that 
coherence is maintained during the calculation. 
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In order to maintain coherence, one must have a small value for e^. At first as one 
decreases e^, the gain in the maximum length number one can factor is roughly inversely 
proportional to the value of e^. However, once e is sufficiently small that the time scale of 
computation for the maximum length which can maintain coherence approaches the inverse 
thermal time 1/T, one reaches a bottleneck. Further reductions in now have little effect 
on the maximum length. The decoherence due to the rapidly increasing time spent in 
computation cancels out the effect of the smaller e^. Thus the thermal time scale 1/T sets 
an effective limit to the time of the calculation, and thus a weak limit on the maximum 
length of the numbers one can compute with. 

If one imagines factoring a 1000 bit number, and one assumes that the quantum factoring 
time can be made to be of order (probably the slowest rate imaginable) , we find that one 
must carry out at least 10® calculation in the thermal time scale. Since the thermal time 
scale for a temperature of IK is of the order of 10~^ sec, this would imply that one would 
have to use a computer which ran at optical frequencies. 

III. OTHER NOISE 

The above coupling to the heat bath is "error free" in the sense that if one is in a 
number eigenstate (ie, is in a state \n >), the system will remain in that state throughout. 
The environment does not cause spin flips. What about the situation in which there is 
also some probability of a state flip- ie of the system making a transition between the two 
eigenstates of cr^? One could approximate this by assuming that the coupling to the heat 
bat is via say 

ae — cos{9)az + sin(^)cr^ 

, with small 6. 

The above analysis is exactly the same for this case, where we replace cr^ everywhere by 
ag. Writing the number eigenstates with respect to ae so that 

\n >e= \nL-i >e ■-■ko >e (31) 
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we have 

\n >= 5^ncos(^)'®"*®'"*sm(^)"*®"^'(-l)^"^®™^)"'|m >e (32) 
The probabihty of remaining in the state \n > under the couphng to the heat bath is then 
Pro6 = 5]^cos(^)2^("''")('>^sin(^)2)(L-,5(n,m)) (33) 

m m 

where S{n, m) is the number of bits in which n and m are the same. Again using the 
arguments above as to the number of terms where the S has a given value, we get 

Proh = {cos{9f + Jsin{eff 

For small 9 this gives 

Proh = (1 - (1 - J)e''f ^ e-^^'-^. (34) 

Thus 6 must be kept very small in order to ensure that the probability of error remains 
small. However we note that the probability of error is vastly suppressed with respect to 
the decoherence probability, which is in accord with the observation that the decoherence 
effects are in general much larger and more rapid than are transition effects. 

This has assumed that the process causing spin flips is the same as the one causing loss of 
coherence in a superposition of the two spin states. In general, the environmental degrees of 
freedom which cause decoherence are not the same as those causing bit flips. I will therefore 
look at the alternative situation in which the single bit Hamiltonian is of the form 

^ ((tti - ei/i(x)a,))2 + {dM^ + (772 - e^h{x)a,)f + (9,02)'(7r3 - e2h{x)ay)f + {dM") 

(35) 

Since we want the single bit decoherence and bit flip probabilities to be small ( or else 
the quantum computer is useless from the start) , I will assume that the are all sufficiently 
small. Furthermore, for simplicity 1 will take £2 = £3, so that the spin flip processes are of 
equal strength. I cannot solve this problem exactly, but since the probabilities are assumed 
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to be very small, one can calculate the transition probability to lowest order in the various 
epsilons. The Hamiltonian can be written as 

H — Ho — ^2 (^J ^iT^i{x)aih{x)dx j h[x)'^dx^ (36) 

where Ui — {az, cr^, ay). The reduced density vector p{t) is given by 

p(t) = Tr (ae'^^'piO) ■ aRre-'''') (37) 

To zeroth order, since H — Hq is independent of a, we have p{t) — p(0). To first order, 
one obtains terms which are linear in the tts and the 0s. However in the thermal state, all 
of these are zero, because the thermal state (of Hq) is symmetric in the fields. To second 
order the results are non-zero. However all of the cross terms €i€j for i ^ j will again be zero 
because the fields are by assumption independent and thus the cross correlations between 
terms linear in each of the fields will again be zero. Thus the only terms surviving will be 
the terms proportional to ej. But each of these terms are independent of the other es. I.e., 
each of these terms are the same as those obtained by setting the other two epsilons to zero. 
These are however just the same as the second order terms calculated above in the first part. 
We thus get 

p{t)i = E,- (l - i Efc {elTr{RT{I{Mt) - Mmdxf))) (38) 
+ i Efc elTriRrifiMt) - (t>m)hdxf)5,k5jk ) P,(0) 

Note that since all of the fields are of the same form and at the same temperature, the 

Tr{RT{J 4>i{t) — (f)i{0))hdx)'^) are the same for all i. 
The probability of bit fiip then becomes 

Probnofiip ~ (1 - < OMt) - 0(O))2|O >)^ (39) 

while the probabihty of decoherence for a state which is the coherent sum over all the integers 
of length L is given by 

Probdecoher ~ (1 - (^(e? + ej) < O|(0(t) - 0(O))2|O > (40) 
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If ei >> 62, the decoherence will again be much more rapid that the probability of "error" 
due to bit flip. 
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IV. CONCLUSION 

Quantum computation places the demand on the system that the coherence of the initial 
state be maintained throughout the computation. In order to maintain this coherence in 
the presence of a heat bath, the reduction in the coupling to the heat bath buys one a pro- 
portional increase in the size of the computation only in the computation can be completed 
within a thermal time scale. For computation times longer than the thermal time scale, a 
decrease in the couphng gives one relatively httle change in the size of the possible coherent 
computation. The thermal time scale thus sets a (weak) hmit on the length of time that a 
quantum calculation can take. 
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